Dynamical pattern formation upon dewetting 
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Dewetting of thin liquid fihns is monitored in situ by 
atomic force microscopy, results are compared with simula- 
tions. The experimental setting is mimicked as close as possi- 
ble using the experimental parameters including the effective 
interface potential. Numerics for the thin film equation are 
based on recently developed schemes which are up to now 
the only methods convergent in all relevant space dimensions. 
Temporal evolution and morphology of experiment and sim- 
ulation are compared quantitatively. Our results explain the 
origin of complex generic patterns that evolve upon dewet- 
ting. 
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For years, experimental and theoretical ||-|T^] 

groups studied spinodal dewetting patterns, yet a quan- 
titative comparison of the results was impossible since 
the potentials of the 'real' experimental system were un- 
known. In a recent study [ p7| , however, we showed how 
to get insight into the underlying potentials by the anal- 
ysis of dewetting patterns. Hence, it is now for the first 
time possible to validate theoretical models also quanti- 
tatively. 

Typically, the evolution of thin liquid films is described 
by highly nonlinear fourth order degenerate parabolic dif- 
ferential equations. In recent years, they became the sub- 
ject both of analytical and numerical studies in mathe- 
matics ||r^-|25[|. It is one peculiarity of these equations 
that they admit globally non-negative solutions. This is 
proved by means of certain integral estimates, called en- 
tropy estimates in the mathematical literature. In order 
to design efficient and convergent numerical tools, one 
is therefore interested in schemes that allow for discrete 
counterparts of the entropy estimate. This has recently 
been accomplished p^-p^ and therefore the method used 
in the present paper can be considered as particularly re- 
liable. 

The stability conditions of liquid films on solid sub- 
strates can be described in terms of the effective interface 
potential W , defined as the excess free energy 

per unit area needed to bring two surfaces from infinity 
to a distance u. If there exists a global minimum of W{u) 
at finite u, the liquid film will not be stable in general; 
rather, it will dewet P-pl|. Dewetting typically proceeds 
by the formation of dry spots, their growth and coales- 



cence, finally leading to a set of droplets on the substrate 
Depending on the sign of W"{u), different rupture 
mechanisms are possible. If W"{u) < 0, the system spon- 
taneously forms dry spots ('spinodal dewetting'), where 
the sites of the holes are correlated reflecting the 

spinodal wavelength Ag. If W"{u) > 0, a nucleus, e.g. a 
dust particle, is needed to induce dewetting. 

In our recent study , we reconstructed the effective 
interface potential W{u) by recording the spinodal wave- 
length as function of film thickness u and by determining 
other system parameters like the contact angle and the 
liquid equilibrium layer. The latter two experimental pa- 
rameters fix the position and the depth of the global min- 
imum of W{u). With the expHcit expression of W{u), we 
now perform new numerical simulations of film rupture. 
In this Letter, we focus on the comparison of experimen- 
tal and numerical results, on the one hand with respect 
to the temporal evolution and on the other hand with re- 
spect to special morphologies of patterns upon spinodal 
dewetting. As we will demonstrate, the appearance of 
so-called 'satellite holes' is not a special pattern limited 
to volatile liquids, as recent work by Kargupta et al. p^ ] 
might suggest, but is a generic phenomenon of dewetting 
liquid films. 

The liquid used in our experiments was atactic 
polystyrene (PS) with a molecular weight of 2 kg/mol 
(purchased from Polymer Labs, Church Stratton, UK, 
M^/M„i — 1.05). This is large enough to assure non- 
volatility, but at the same time small enough for the poly- 
mer melt to be approximately a Newtonian fluid. This 
is essential for comparison with simulations, in which 
viscoelastic effects were neglected. Films were prepared 
by spin casting a toluene solution (Selectipur toluene, 
Merck, Germany) of PS onto polished, oxidized Si wafers 
(Silchem GmbH, Freiberg, Germany). Before coating, 
the wafers were thoroughly cleaned using standard pro- 
cedures HJl^,^. The oxide layer and PS film thickness 
were determined by ellipsometry (Optrel GdBR, Berlin, 
Germany). The surface of the wafers, called 'SiO-wafers' 
below, consists of amorphous, 191(1) nm thick silicon 
oxide On such a substrate, PS films are unsta- 

ble up to a film thickness of about 300 nm ||2^. The 
symmetry-breaking mechanism of dewetting was moni- 
tored by atomic force microscopy (AFM) (Nanoscope HI, 
Digital Instruments, Santa Barbara) in tapping mode"^*^. 
Annealing took place on top of the AFM sample holder 
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which enabled us to follow the dewetting process at a 
real time scale. We therefore continuously scanned the 
sample and recorded one scan in every 60 s. Scanning pa- 
rameters were carefully adjusted not to affect the liquid 
layer |§. 

Figure 1 depicts on the left side a series of AFM cross 
sections of a 3.9(1) nm thick PS film dewetting a SiO- 
wafer at an annealing temperature of Ta — 53 °C. The 
film is unstable and holes are generated due to spinodal 
dewetting of the film |jl^,^. The characteristic wave- 
length Xs can hardly be observed in the cross sections, 
since the scan size is too small. Before annealing, the sur- 
face of the PS film is smooth, exhibiting an rms rough- 
ness of below 0.3(1) nm. As revealed from the online 
AFM cross sections, first holes are detected after roughly 
1100 s. Looking at the entire (1.5 /im)^ AFM scan, 
though, holes pop up within a time interval of about 
750 s < t < 4000 s, due to slight thickness variations 
of the polymer film. Hence, the cross sections shown in 
Fig. la are a kind of 'random sample' for the spinodal 
break-up of a hole. On the bottom of the holes, a 1.3 nm 
thin layer of PS is left, as was measured by X-ray reflec- 
tometry. It is identified with the equilibrium thickness of 
PS on SiO 0. During the entire dewetting process, the 
contact angle of 7.5(5)° was found to remain constant. 
Due to mass conservation, the material removed from the 
inner side of the hole is accumulated at the perimeter of 
the hole. 

The effective interface potential for that system can be 
written as IITtII 



where the stabilizing term W+ and the destabilizing term 



Jir/ \ ^ ^SiO 



(1) 



where e denotes the strength of the short range part of 
the potential, e = 6.3(1) • 10"''^ Jm^, and Asio is the 
Hamaker constant of PS on SiO, As^o = 2.2(4) • 10"^° J. 
The global minimum in W{u) is at u=1.3(l) nm reflecting 
the equilibrium thickness of a PS film on top of SiO- 
wafer. W{u) enters the equation of motion of the liquid 
film. 



du 

rj— — div(m(u)Vp) = 0, 

via the pressure p which is given by 
p = -a-Au + W'{u), 



(2) 



(3) 



where a and rj denote the surface tension and the viscos- 
ity of the film, respectively. A no-slip condition at the 
solid-liquid interface entails m{u) — ^u'^- 

For numerical purposes it is important to distinguish 
between stabilizing and destabilizing terms. Therefore, 
we decompose 



W- are given by eu and 



respectively. 



In general, smooth solutions of system (|^)-(^) are not 
expected to exist. Instead, we are looking for so called 
weak solutions u,p that satisfy 

nr]Ut(j)+ [ ( m{u)\/pV<l) = Q (5) 

[ ( p^Ij^ ( [ crVuV^+ [ [ W'{u)^/j (6) 
Jq Jn Jo Jn Jo Jn 

for all V, e C^{iO,T) x fl). Note that smooth solu- 
tions are also weak solutions. Moreover, (^)-(^) are the 
starting point to develop finite element schemes (cf. |^ ) . 

For simplicity, we sketch the main idea in one spa- 
tial dimension: Let the space interval fl be divided in 
subintervals [cc^-i, Xi], 1 < i < D of equal length h. Sim- 
ilarly, assume a time discretization with nodal points tk , 
< A: < iV to be given. We look for discrete solutions 
in the linear finite element space V^. A canonical ba- 
sis of V'^ is given by functions tpi which are linear on 
each subinterval and which satisfy ipi{xj) = Sij. Our 
discretization is as follows: 

For given g V'^ , find a sequence {U^,P'^) £ y'' x 
V^,k^O,...,N -I such that 

IL{ljk+i _ [/fc^ $),^ + {M{U^+^)d^P^+^,d^^) = (7) 

Tk 

{p''+\-^)h=(7{d,u''+\d.,^) 

+ {W'^{U''+^) + W'_{U''),^)h 

for all e F'*. 

Here, (., .) denotes the scalar product in the Lebesgue- 
space L^(r2) and (., .)h is the scalar product of lumped 
masses, which is an approximation of (., .) given by 
{ipi,<Pi)h = Z]j('/'i>Vi) and {ipi,Lpj)h \ii ^ j. Cru- 
cial is the choice of the discrete version M of the mobility 
TO. The right strategy is to define M on each interval 
a^i+i] as 



M{u) = \ ^ J if ^(^') ^ 

m{U{xi)) otherwise 



(9) 



where denotes the mean value integral. 

The procedure discussed above naturally endows us 
with discrete versions of well-known integral estimates 
which read 



ry / {a\Vu{T)\^ + W{u{T))) 



TO(li)|Vp|^ 



Jn 



W{u) = W+{u) + W^{u), 



(4) 



<77 / {a\Vuo\^ + W{uo)) (10) 
Jn 
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and 

7j f G{u{T)) + a r[ \Au\^ + W:i{u)\\/u\^ 
Jn Jo Jn 

<V [ G{uo) + K n \Wu\\ (11) 
Jn Jo Jn 

where G denotes a second antiderivative of m^^. They 
are obtained from by setting (jj ~ p and (t>{x, t) — 
G'{u{x,t)), respectively. 

Combining these estimates, it can be shown that dis- 
crete solutions are positive and that in addition both 
U and P converge strongly with respect to the topol- 
ogy of the Lebesgue-Sobolev-space L'^{{0,T), H^^'^{n)). 
Moreover, U converges strongly to u in the Holder-space 

The right-hand side of Fig. 1 depicts numerical results 
of thin film dewetting. The sequence makes use of the 
effective interface potential as reconstructed for PS film 
on top of a SiO-wafer shown in Eq. (1), the viscosity 77 = 
12000 Pa s [|0) and the surface tension a= 30.8 mN/m. 
In the numerical simulation we used 1000 grid points and 
an initial film of height 3.9 nm with a slightly roughened 
surface. 

The evolution of a hole in the simulation resembles 
the experimental AFM scans. The absolute time scale is 
in the same order of magnitude and the morphological 
development of a hole in experiment and simulation is 
very similar. Deviations can be observed especially in the 
beginning, where the initial conditions of the film surface 
play a significant role. In the experiments, the rupture 
time varies in a time interval of 750 s < t < 4000 s which 
is comparable to the time scale in the simulation. 

If we start the simulation, however, with an 'artificial 
step' between the equilibrium layer and the original film 
thickness of 4.9 nm, which is now slightly thicker than 
before, it generates a pattern like the one shown on the 
right-hand side of Fig. 2. The step transforms rapidly 
into a heap that relaxes via an undulation into the origi- 
nal film thickness. With time, the dips of the undulation 
grow deeper and the first dip reaches the substrate after 
about 400 s. As the other dips follow to break through 
to the substrate, a series of consecutive holes appears, 
which are generated by a kind of 'cascade effect'. 

A strikingly similar morphology can be observed ex- 
perimentally, too. Fig. 3a-d depict a temporal series 
of (10 nm)2 AFM scans of a 4.9 nm thick PS(2k) film 
on a SiO-substrate. First, single holes appear that soon 
grow in size. From a certain size onwards, they are sur- 
rounded by a second row of holes, followed later by a third 
one and so on. We term this pattern 'satellite holes'. 
As compared to the experimental situation reported be- 
fore, the film now is 1 nm thicker. The film can still 
dewet spinodally, yet the growth time of the amplitude 
of the preferred mode is about a factor of three longer 
|p|,p3, therefore dewetting by nucleation is the quicker 



way of dewetting. The onset of the spinodal wavelength 
yet can be seen by the wavy pattern on the film matrix, 
cf. Fig. 3b and c. 

A comparison of simulation and experiment is shown 
in Fig. 2. The experimental profiles depicted in the left 
side of Fig. 2 are horizontal cross sections of the AFM 
scans of Fig. 3, taken at the level indicated by the arrow; 
t = is defined as the time when first small holes can 
be detected in the cross section. Depending on the direc- 
tion the cross sections are cutting through the AFM scan, 
first satellite holes can be seen in a time interval between 
300 s and 1600 s. The time scale of experiment and sim- 
ulation hence is in the same order of magnitude. Slight 
differences are to be expected, since i) the viscosity of the 
experimental system is known only by roughly a factor 
of two and ii) effects due to slippage of the polymer melt 
arc not taken into consideration in the simulation. Yet 
in both films, an undulatory behavior of the rim profile 
can be observed leading to satellite holes. This is clear 
evidence that this dewetting pattern does not require a 
volatile liquid, but is generic to a spinodally dewetting 
system. 

In a recent study, we showed that an undulatory rim 
profile can be observed on polymer films of thicknesses 
up to at least 60 nm, if the polymer melt is below the 
entanglement length . There we showed that theoret- 
ically a disturbance in a Newtonian liquid should always 
decay via an undulation and that only the viscoelastic 
properties of the liquid may stabilize the system. The 
observation of an undulatory rim profile in the simula- 
tion, where no viscoelastic properties were considered, is 
in accordance with this assertion. 

A further interesting observation in both experiments 
and simulations is that a dip of an undulation only leads 
to a hole if the initial film thickness is below about 7 nm. 
The reasons for this are speculative at the moment: It is 
possible that the driving force for a dip to grow deeper 
falls below a certain value for films thicker than 7 nm. It 
is also possible that there is an interplay of the dewet- 
ting velocity of the three phase contact line, as set by 
the global minimum of W(u), and the viscous flow of 
the melt inside the dewetting front. Further experiments 
and new simulations will surely shed more light into this 
discussion. 

To sum up, we may say that lubrication approximation 
is an appropriate method to model the dewetting of liquid 
films also beyond the time of the first hole formation. We 
thus have a powerful tool to test the influence of single, 
yet different parameters onto the dewetting pattern and 
its temporal evolution. A successful example was given in 
the form of the 'cascade' effect leading to satellite holes. 
It was shown that this is a generic pattern of a spinodally 
dewetting fllm. We intend to perform numerical studies 
of the three dimensional system in a forthcoming paper. 
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FIG. 1. Cross sections of a spinodal rupture scenario in a 
3.9 nm thick PS(2k) film on a SiO-substrate as seen in situ 
by AFM (left) in comparison with a simulated one (right). 

FIG. 2. 'Cascade effect' leading to satellite holes in a 
4.9 nm PS (2k) film on a SiO-substrate; AFM cross sections 
(right) and numerical simulation for rj — 1200 Pa s (left). 

FIG. 3. Pattern formation in a 4.9 nm PS(2k) film on a 
SiO substrate exhibiting 'satellite holes' a) t = s, b) 866 s, 
c) 1303 s, and d) 3060 s. Horizontal cross sections of the scans 
are taken at the height of the arrow and are shown at the left 
side of Fig. 2. 
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